Bowen-York trumpet data and black-hole simulations 
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The most popular method to construct initial data for black-hole-binary simulations is the punc- 
ture method, in which compactified wormholes are given linear and angular momentum via the 
Bowen-York extrinsic curvature. When these data are evolved, they quickly approach a "trumpet" 
topology, suggesting that it would be preferable to use data that are in trumpet form from the 
outset. To achieve this, we extend the puncture method to allow the construction of Bowen-York 
trumpets, including an outline of an existence and uniqueness proof of the solutions. We construct 
boosted, spinning and binary Bowen-York puncture trumpets using a single-domain pseudospectral 
elliptic solver, and evolve the binary data and compare with standard wormhole-data results. We 
also show that for boosted trumpets the black-hole mass can be prescribed a priori, without recourse 
to the iterative procedure that is necessary for wormhole data. 

PACS numbers: 04.20.Ex, 04.25. Dm, 04.30.Db, 95.30.Sf 



I. INTRODUCTION 

Numerical solutions of the full Einstein equations for 
the last orbits and merger of compact binary systems are 
important for the developing field of gravitational-wave 
astronomy. In the case of black-hole binaries, long-term 
simulations became possible in 2005 0, 0, Q , and within 
the last few years the field has developed to the point 
where the gravitational-wave (GW) signal from such sys- 
tems can be calculated to essentially the required accu- 
racy of current GW detectors and work is underway 
to incorporate these results into GW searches Q. How- 
ever, only a small fraction of the full black-hole-binary 
parameter space has yet been studied @, and its full 
exploration will require yet more accurate and efficient 
numerical simulations. The first step in any simulation 
is the production of initial data, and these determine in 
part the accuracy and physical fidelity of the final simu- 
lation; that is the focus of this paper. 

The 3+1 approach to solving Einstein's equations con- 
sists of specifying initial data (the metric and its time 
derivative on one constant-time slice of spacctimc), and 
then evolving that data forward in time. Valid initial 
data satisfy a set of constraint equations, and a given so- 
lution to the constraints will represent a certain physical 
situation in a certain set of coordinates. We are then 
faced with the problem of finding constraint-satisfying 
data that both represent the physical situation we wish 
to simulate (in our case two black holes following non- 
eccentric inspiral) and are in a suitable set of coordinates. 

The most widely used method to evolve black-hole- 
binary initial data is the moving-puncture method @, Q , 
which involves a modification of the Baumgarte-Shapiro- 
Shibata-Nakamura (BSSN) @, § formulation of the 3+1 
ADM- York Einstein equations [9, [l(| combined with the 
"1+log" [ll| and "f -driver" gauge conditions [13, El- 

As the name suggests, the data that are usually evolved 
with this method are puncture data [HI, whereby black 



holes are represented on the numerical grid by compact- 
ified wormholes. However, when these data arc evolved 
using the standard moving-puncture method the numer- 
ical slices lose contact with the extra asymptotically flat 
wormhole ends, and quickly asymptote to cylinders of fi- 
nite areal radius located within the horizon of each black 
hole. That the data evolve to these "trumpets" was re- 
alized in [l5j . in which an analytic stationary trumpet 
endstate was derived and shown to agree with numerical 
results. 

That work suggested a new form of initial data, based 
on trumpets. It was shown in [l6| that maximally- 
sliced trumpet data can easily be constructed numeri- 
cally based on the solution first presented in [ijj, and 
that these data are indeed time independent in a moving- 
puncture simulation. These data represent the first non- 
trivial test solution for most current black-hole evolution 
codes. It was later shown in [l~8[ that an implicit form of 
the same solution could be constructed analytically, and 
in [lj| the corresponding solution for the 1+log-sliced 
case was found. 

We have presented a detailed study of Schwarzschild 
wormholes and trumpets in [lj|, with a focus on con- 
structing and evolving Schwarzschild trumpet puncture 
data. This work extends that study to boosted, spinning 
and binary tr ump ets. As we described in the conclud- 
ing section of [l9(, ideal binary puncture data will be in 
trumpet form, 1+log-sliced (or satisfy whatever slicing 
condition is ultimately used to evolve them), and repre- 
sent true boosted Schwarzschild or Kerr black holes (i.e., 
will be free of the junk radiation that plagues all current 
binary simulations). As a first step in a larger research 
programme to attempt to achieve that goal, we deal here 
with only the first point in our list of requirements: that 
the data be in trumpet form. The data we construct 
will not meet any of the other requirements: they will be 
maximally (not 1+log) sliced, and they will be confor- 
mally flat, meaning that they include essentially the same 
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junk radiation as standard puncture data. As such, this 
work is a proof-of-principle exercise that demonstrates 
that is feasible to produce binary trumpet data. Along 
the way a number of new issues arise that are not present 
in the wormhole case, and variants of these issues may 
recur in efforts to produce yet more general data. 

We will start with a brief summary of wormholes, 
trumpets and punctures in Sec. [TH then discuss in some 
detail the maximal slicing case in spherical symmetry in 
Sec. IIIII as an example for our analytical setup to con- 
struct trumpet data and for our numerical algorithm. 
We then extend the trumpet-puncture construction to 
boosted (Section HVj) and spinning (Section |V| Bowen- 
York black holes, and provide an outline of a proof for 
both existence and uniqueness of these solutions. In Sec- 
tion I VII we estimate the junk-radiation content of these 
initial-data sets, before moving on to binary data in Sec- 
tion |VlTl The ultimate goal is of course to produce data 
that can be used in black-hole-binary simulations, and 
in Section [Villi we evolve a binary data set and compare 
with the corresponding standard wormhole-puncture re- 
sults. We close with a discussion on the next steps to 
producing optimal initial data for moving-puncture sim- 
ulations. 



II. BACKGROUND: A BRIEF SUMMARY OF 
WORMHOLES, TRUMPETS AND PUNCTURES 

A. Wormhole puncture data 

Consider a constant-time slice of the Schwarzschild 
spacetime. Write the standard Schwarzschild solution 
in isotropic coordinates, i.e., 



ds 2 




dt 2 +ip 4 {dr 2 + r 2 dfl 2 ), (1) 



and the isotropic coordinate r is related to the 
Schwarzschild areal radial coordinate R by 



R 



ip 2 r, 



(2) 



and ip = 1 + M/2r is a conformal factor. Now the 
data (7y , Ky ) on any t = constant slice are given by 
jij = 4) A r)ij (where ijij is the flat-space metric in the cho- 
sen coordinate system) and = 0. The fact that the 
physical spatial metric can be related to the flat-space 
metric using only the conformal factor indicates that the 
solution is conformally flat. 

We see immediately from Eqn. ([2]) that the slice does 
not reach the physical singularity at R = 0, or even pen- 
etrate the black-hole horizon at R = 2M. In fact, the 
coordinate range r G [0, oo] contains two copies of the 
Schwarzschild spacetime exterior to R = 2AI: one copy 
in r G [0, M/2] and the other in r G [M/2, oo]. These co- 
ordinates therefore represent the exterior Schwarzschild 
spacetime as a wormhole, and this is most clear when 



viewed in an embedding diagram like that shown in Fig. 1 

of [3. 

The advantage of these slices for numerical relativity 
is that the entire exterior space can be represented on 
R 3 without any need to deal explicitly with the physical 
singularity of the black hole, or to "excise" any region 
of the computational grid. The point r = 0, which is 
commonly referred to as a "puncture" (l4| . represents a 
second copy of spatial infinity, but the solution is well- 
behaved there, except for the conformal factor ip, which 
diverges as 1/r. 

We can write initial data for multiple Schwarzschild 
black holes simply by modifying the conformal factor to 
ip = mj/(2rj), where the m, parametrize the mass 

of the ith black hole, and the ith puncture is located at 
Vi = po| . Furthermore, one may imbue these black 
holes with linear and angular momentum by providing a 
non-zero extrinsic curvature. If we retain the property 
of conformal flatness and choose the extrinsic curvature 
to be trace- free (K = 0, or maximal slicing), then there 
exist solutions of the momentum constraint for boosted 
and/or spinning black holes; these are the Bowcn-York 
solutions 21] . The solution is provided only in the con- 
formal space, and is related to the physical extrinsic cur- 
vature by 



R~ij — V' ^3 > 



(3) 



where here Aij is the Bowen-York solution. Now, how- 
ever, the conformal factor is not known analytically, and 
can only be found by solving numerically the Hamilto- 
nian constraint, 



~9 1 

V 2 tp + -ip 



' 7 A i3 A tj = 0. 



(4) 



The most convenient way to solve (U) is by the "puncture 
method" [l4| , which is to realize that the solution can be 
constructed conveniently in terms of a (typically small) 
correction u to the Brill-Lindquist solution, 



4> = i + 



. 2^ + ' 



(5) 



Since the Brill-Lindquist conformal factor is in the kernel 
of the flat-space Laplacian, the Hamiltonian constraint is 
now an equation for the correction function u: 



V 2 u + -ip~ 



0. 



(6) 



Furthermore, the function u is sufficiently regular over all 
of i? 3 that © is in the form of a nonlinear elliptic equa- 
tion that is straightforward to solve by a number of stan- 
dard methods. This approach is used to construct the 
majority of black-hole-binary initial data used in current 
numerical simulations, and the elliptic solve is performed 
either with mesh-refinement finite-difference solvers [22| 
or, in most cases, by an elegant single-domain spectral 
approach [23| , which we will adopt for the work presented 
here. 
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A key property of the Bowen-York family of solutions 
to the momentum constraint is that the values of the mo- 
mentum and angular momentum of the spacetime (and 
thus in some sense the momenta and spins of the black 
holes) can be prescribed before solving for the conformal 
factor. 

Two further properties of these data deserve particular 
attention here, and we will return to them when compar- 
ing these wormhole puncture data with our new trumpet 
puncture data in the following sections. 

First, Bowen-York black holes are conformally flat, 
which is not the case for either a true boosted 
Schwarzschild black hole, or a Kerr black hole, or a 
boosted Kerr black hole. Since it is usually a boosted 
Schwarzschild or Kerr black hole that we really wish to 
describe, these data are often described as the desired 
physical objects plus some "junk"; the junk represents a 
valid part of a solution of Einstein's equations, but it is 
not a part that we would expect to occur physically, and 
can be interpreted as unphysical gravitational wave con- 
tent. As the data evolve forward in time, the junk either 
falls into the black hole or radiates away, quickly leaving 
precisely the physical situation that was intended in the 
first place, albeit with slightly different physical param- 
eters. 

In practice (i.e., in black-hole-binary simulations) this 
junk radiation causes two problems. One is that it in- 
troduces noise into the numerical simulation, which can 
affect the numerical accuracy. This point is clearly illus- 
trated in [24[. The other is that it limits the physical 
black-hole spin that can be achieved. When the spin 
angular momentum of the Bowen-York black hole is ex- 
tremely high, most of the angular momentum manifests 
itself as junk, and after that junk has either fallen into the 
black hole or radiated away, we are left with a Kerr black 
hole that has spin no higher than a/M = S/RI 2 < 0.93 
l26l . H3] ; we will confirm this with high- precision nu- 
merical simulations, bounding the final Kerr parameter 
at a/m < 0.929. This property of Bowen-York data 
preclude their use to study very highly spinning black 
holes, which may in fact be the most common astrophys- 
ically [H,[2^,[30|, and we must turn to other types of data 
— see, for example, [3l| for the construction and evolu- 
tion of spinning but non-boosted puncture data, and [27| 
for non-conformally-flat black-hole initial data where the 
interior of the black hole is excised. 

The other property of Bowen-York puncture data that 
we want to highlight is the calculation of the black-hole 
mass. Having produced data for two black holes, we 
would like to know what their masses are; although the 
parameters m, parametrize the black-hole masses, the 
black-hole mass equals the mass parameter only in the 
case of a single Schwarzschild black hole, i.e., the origi- 
nal Schwarzschild solution in isotropic coordinates. 

In any other case, we typically estimate the black-hole 
mass by two methods. One is to calculate it from the 
area of the apparent horizon. This requires that we first 
locate the apparent horizon, which can be computation- 



ally expensive (although fast and efficient solvers exist, 
for example [32j|). The other method is to make an inver- 
sion transformation at each puncture and calculate the 
ADM mass at that black hole's extra asymptotically flat 
end, and to treat this quantity as the black-hole mass. 
For a binary system, this mass estimate is given by 

Mi = rrn \l + uo,< + -^1) , (7) 

where Dij is the coordinate separation between the two 
punctures, and uo.i is the value of the correction func- 
tion u at the ith puncture. Remarkably, this expression 
is found to agree within numerical error with the mass 
calculated from the apparent horizon [33|, [H| , although 
we will see in Section IVT1 that this can only be expected 
to hold for boosted black holes, or black holes with small 
spins. 



B. Trumpet puncture data 

Bowen-York puncture data were first constructed long 
before stable numerical simulations of black-hole bi- 
naries were possible, and were useful in both mathe- 
matical relativi ty [35l. l36l . [3tJ and in studies of initial 
data (Til. l38l. [39l. |40|. l4l| . However, with the advent of the 
moving-puncture method 0, H[ it was found that worm- 
holes may not be the most suitable topology for black- 
hole initial data. 

In a moving-puncture simulation, the numerical slices 
quickly lose contact with the extra asymptotically flat 
ends, and instead asymptote to cylinders of finite areal 
radius [H, M, EI 1MB, or "trumpets" . This sug- 
gests that it would be more natural to construct initial 
data in trumpet form from the outset. 

To date this has only been done for a single 
Schwarzschild black hole. The question addressed in this 
paper is, How can we generalize the wormhole puncture 
procedure to produce trumpet punctures for black-hole 
binaries? For a single maximally sliced Schwarzschild 
black hole, the trumpet data can be put in a form simi- 
lar to the wormhole isotropic coordinates, where now the 
conformal factor behaves as ip ~ y / 3M/2r near the punc- 
ture. However, the full conformal factor is not known an- 
alytically (except as an implicit equation in terms of the 
Schwarzschild radial coordinate R) [TU, [l9[ . This means 
that it is not straightforward to superpose two trumpets 
as with the Brill-Lindquist solution in the wormhole case. 
And it is not obvious how the introduction of the Bowen- 
York extrinsic curvature (which, if we retain conformal 
flatness and maximal slicing, remains a valid solution of 
the momentum constraint), affects the behavior of the 
conformal factor near the puncture, or the physical prop- 
erties of the data. Finally, without the presence of extra 
asymptotically flat ends, we lose the simple procedure to 
estimate the black hole's mass from Eqn. ([7]). These are 
the issues that we address in this work. 
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In Section Hill we describe in more detail the maximal 
Schwarzschild trumpet, and use it to illustrate our more 
general method for producing single-trumpet data. 



In order to solve the Hamiltonian constraint, we start 
with an ansatz for ip that includes the required asymp- 
totic behavior. We write the full conformal factor that 
solves the Hamiltonian constraint as 



III. MAXIMAL SCHWARZSCHILD TRUMPET 

A. Constructing a conformal- factor ansatz for 
trumpet data 

The basis of this work are data that represent a max- 
imal slice of the Schwarzschild spacetime with a trum- 
pet topology. The first hints of this representation of 
Schwarzschild were given by Estabrook et al. [171 ] in 



1973, but it wasn't until the development of the moving- 
puncture method 0, Q in 2005, and a subsequent un- 
derstanding of the dynamical behavior of the numeri- 
cal slices [ 1 5| in that method, that it was realized that 
the maximal Schwarzschild trumpet could be expressed 
in a simple form [T(jj |. and could in turn be written in 
the "puncture" isotropic coordinates suited to moving- 
puncture simulations [lfjl . [l8j . 

For a single Schwarzschild black hole with mass M, the 
conformal initial data in Cartesian coordinates are 



K 



a 



ft' 



— (3riinj - Sij) , 
0, 



1 - 



2M C 2 



x l aC 



where C = -y/27/16M 2 , R is the Schwarzschild radial 
coordinate, r — (x 2 + y 2 + z 2 ) 1 / 2 is the isotropic radial 
coordinate, and n% = Xi/r is the outward-pointing nor- 
mal vector. All that remains to fully specify the initial 
data is a valid conformal factor tp that maps these data 
to the physical space, i.e., 



Hi 



ip%j 



1)- 
ip 2 r. 



The conformal factor must satisfy the Hamiltonian con- 
straint and asymptote to ^ — > 1 as r — > oo. A numerical 
solution of the Hamiltonian constraint for these data was 
first presented in (lfj| , and an analytic solution (albeit an 
implicit solution in terms of i?, not r) given in (l8j . 

To illustrate the method that we will use for more 
general cases, and to test our elliptic solver, we will 
again solve the Hamiltonian constraint numerically. Our 
boundary conditions are that ip —* 1 as r — > oo, and 
tp ~ v/3M/2r as r — > 0; the latter condition ensures that 
we have a trumpet topology. 



Ip = 1p s + U, 

where ip s incorporates the desired asymptotics. 
Hamiltonian constraint for this problem is 



\7 2 u 



7 A - ■ A l 3 



(8) 



The 



(9) 



where V 2 represents the Laplacian with respect to the 
flat background metric, and it is understood that A^j — 
Afj, although this is the form of the Hamiltonian con- 
straint that we will deal with for all choices of Aij 
throughout this paper. 

One easy way to incorporate the asymptotic behavior 
is to apply weight functions to the two asymptotic con- 
ditions, 



3AI 



ip s (r) = wi(r)\j + w 2 {r) 



such that 



u>i(0) = 1, w\{oo) = 0, 
w 2 (0) = 0, w 2 (oo) = 1. 

The weight functions we choose are 
w 2 (r) 



1 + r 4 

^4 



1 + r 4 

These have the property that at each end of the slice 
the conformal factor's lowest-order deviation from the 
required behavior is at fourth order. 

Consider now the behavior of the conformal factor near 
the puncture. We assume the leading order terms to be 
of the form 



A 

,1/2 



Br r - 



(10) 



If we insert this ansatz into the Hamiltonian constraint, 
we have 



-r 1 



A 



4 r 5/; 



+ Bn(n + l)r r ' 



81 M 4 
64r 6 
81Af 4 
64A 7 r 5 / 2 X 

1 _ 7 JL r n+i/2 
A 



where we have expanded about r = on the right-hand 
side. Equating coefficients of r, we find that A = ^/3M/2 
(as we expect). We also find that for a consistent solution 
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n = \/2 — 1/2 = 0.9142... and B remains undetermined. 
We therefore see that divergent terms near the puncture 
do exactly cancel, and the next-to-leading order term 
goes to zero. However, this next-to-leading order term 
goes to zero with a non-rational power of r (which was 
also noted in |45j). and this may limit the accuracy of a 
spectral solution to ©■ If this is the case, we may also 
include the yv2-i/2 behavior into our ansatz. 

An implicit solution of tp in terms of the Schwarzschild 
radial coordinate R is given in [l8T ]. as is an implicit 
solution of r(R). If we combine these as (ip(R) — 
v / 3M/2r(R))/r(R)^ 2 ~~ 1 / 2 , and take the limit as R -> 
3M/2, we can determine the coefficient B in our ansatz 
above. We find that 



B 



3M\ 



3/2 



M 



3MV 
2^2/ 



-1—V2 



fll) 



If necessary, we may now use 



Mr) = Mr) (A + Br^- 1 ^ + w 2 (r), (12) 

as the ansatz in our numerical solution of the Hamilto- 
nian constraint. 

To summarize, we have two choices of conformal factor 
ansatz that we may adopt, and which we denote by, 



ip = ipsi (r) + u 
ip = ips2(r) + u 



3M 



w i( r )\l — + w 2(r) +u, (13) 



1 (r)( A +Br V2-i/2 



+w 2 (r) + u . 



(14) 



B. Numerical solution of the Hamiltonian 
constraint 

In order to solve the equations numerically, we have 
written a code to solve systems of nonlinear elliptic equa- 
tions with general finite difference methods in three spa- 
tial dimensions. In this work we will only utilize this 
solver with pseudospectral discretizations, representing 
the solution by Fourier series in (periodic) angular coor- 
dinates, and as Chebyshev polynomials otherwise. The 
solver has been developed as a Mathematica package, 
it uses the Mathematica LinearSolve function with a 
Krylov method and ILU prcconditioner to solve Linear 
systems, and Newton iteration to deal with nonlinear- 
ities. This approach has allowed us to develop a very 
flexible spectral elliptic solver from scratch, in order to 
achieve good performance even for the larger grids we 
use in this paper. We consistently use sparse matrix ob- 
jects and generate compiled code using Mathematica's 
CompiledFunction for certain key functions which oper- 
ate on individual matrix elements. 

The elliptic solver uses compactificd coordinates 
(X, Y, 0), with X G [-1,1], Fe (-1,1) and G (-< 7r,7r). 




FIG. 1: The L2 norm of the error in the solution function u for 
a maximal Schwarzschild trumpet. The dashed line shows the 
error when using the ansatz (|13[) . while the solid line shows 
the error when using the ansatz ()14l) , which includes the next- 
to-leading order behavior in the conformal factor near the 
puncture. 



In all cases that involve a single black hole, we trans- 
form to these coordinates from spherical polar coordi- 
nates with r = (1 — X)/(l + X) and Y = cos(0), so that 
X = — 1 corresponds to r — > 00 and X = 1 corresponds 
to r = 0. In order for the coefficients of the Laplacian op- 
erator to be sufficiently smooth over the entire domain, 
the entire equation is weighted by a factor 



MX,Y,< 



(1+X) 3 (1-Y 2 ) 



(15) 



The accuracy of the numerical method is demonstrated 
in Fig. [TJ which shows the L 2 norm of the error between 
the numerical and analytic solutions as a function of the 
number of collocation points N. (The same number of 
points is chosen in each direction, although since this 
solution is spherically symmetric, the solution varies only 
along the X direction.) It is clear from Fig. [T]that the 
spectral convergence is lost for N > 20 when the ansatz 
ipal is used, but remains up to at least N = 48 where the 
next-to- leading order behavior is included in i\) s2 . 

The numerical solution u is shown in Fig. O Solutions 
using both the ip s i and i\) s2 ansatze are shown. The sec- 
ond panel in the figure zooms into the region near the 
puncture. In this figure the solution was produced us- 
ing the ansatz with ip S 2- We can see that the function 
smoothly approaches zero at the puncture, and is well 
resolved by the numerical method. The i/^-based so- 
lution is not well resolved near the puncture and is not 
included in the second panel. The data in this plot are 
from solutions with TV — 52 collocation points. 



IV. SINGLE BOOSTED BOWEN-YORK 
TRUMPET 

We now consider a single trumpet with linear momen- 
tum. 
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FIG. 2: The correction function u for the maximal 
Schwarzschild trumpet, for the solution with N = 52 colloca- 
tion points. The solution using the ip s i ansatz (|13|) is shown 
with a dashed line, and the solution using the ip S 2 ansatz (|14|) 
is shown with a solid line. The second panel zooms into the 
region near the puncture, to illustrate that the ^,,2-based so- 
lution smoothly approaches zero there, and is well resolved; 
the i/> s i solution is poorly resolved in this region. 



To do this we add to the conformal extrinsic curva- 
ture the Bowen-York solution for a single black hole with 
linear momentum P l , 

3 

A% Y (r; P) = ^ (Pirij + P 3 n % - (% - mn j )P k n k ) , 



so that the total conformal extrinsic curvature is 



(16) 



(17) 



The superscript "BYT" is a reminder that this is a 
Bowen-York trumpet. 

Note the asymptotic behavior of A® Y as r — > 0: it 
diverges as 1/r 2 . Since the trumpet extrinsic curvature 
diverges faster, as 1/r 3 , it dominates the Hamiltonian 
constraint near the puncture, and so determines the be- 
havior of the solution. In particular, this means that the 
trumpet form ip ~ ^/3M/2r remains. 

Consider the general form of AijA 1 ^ for the boosted 
case (in the following we will suppress the angular de- 
pendence of the functions for simplicity): the contribu- 
tion from the Schwarzschild trumpet diverges as r -6 , the 
Bowen-York contribution diverges as r -4 , and the cross 



terms diverge as r , so we have 

J rp^± ipO y>0 



(18) 



Near the puncture, we can write the inverse conformal 
factor term as 

t/T 7 = (Ar- 1/2 +u)- 7 

= A- J r 7 / 2 (l + ur^ 2 /A)- 7 

= A- 7 r 7 / 2 -7A- 8 ur A + 0{r 9 ' 2 ), 

where A is the same quantity that was introduced in 
Eqn. (|10p . We can now write out the source term of the 
Hamiltonian constraint as 



< A - ■ A l 3 



E 

i=0 



D, 



r l/2+i 



E 



uD' 



(19) 



The D2 term is the one that diverges as r~ 5 / 2 and is 
canceled by a corresponding term from the Laplacian of 
ip s , as described in Section IlIII The remaining terms all 
result in contributions to u with positive powers of r, and 
which therefore go to zero at the puncture, except for the 
D' 2 term, which can in principle lead to a contribution 
that diverges as lnr. We note that such a term also 
appears in the Schwarzschild case (with our choice of 
ansatz), but there we know that u = at the puncture, 
and so none of the D[ terms contribute to the solution. 
Fortunately, we will see in the existence proof that we 
present below that the same is true in the boosted case. 
In the coordinates of our elliptic solver, the puncture r = 
is located on the entire coordinate plane X = — 1, and 
so there we can simply impose that either u = or u = 0, 
and thus prevent the solver from producing unphysical 
divergent terms. 

Before proceeding, we will show that solutions to this 
problem exist and are unique. Note that while construc- 
tion of a numerical solution gives evidence for the exis- 
tence of a solution to the continuum equations, unique- 
ness is not easy to verify numerically, and an analytical 
proof is highly desirable. While the uniqueness proof is 
general, the existence proof requires a more detailed anal- 
ysis of the Hamiltonian constraint, and in the spinning 
case we will deal with only a single trumpet (we do how- 
ever expect that the same procedure can be generalized 
to multiple spinning and boosted black holes). 

We first prove uniqueness. Assume we have two posi- 
tive solutions, tpi and Subtract the equations to get 



v 2 (Vi-^ 2 ) 



l -A l3 A^(^ 7 



We assume V12 = ipi — i>2 goes to zero at both ends 
(they satisfy the same boundary conditions, and we saw 
in the preceeding discuss that there are no other diver- 
gent terms in the solution). If ^2 is not identically zero, 
it must have a positive maximum or a negative minimum. 
Neither of these is compatible with the equation (leading 
to different signs on the left and right hand side). 
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We now provide the outline of an existence proof. The 
maximum principle tells us that a solution, if it exists, 
cannot have an interior minimum. As r — > oo our bound- 
ary condition is that iji — ► 1, and so a solution, if it exists, 
satisfies ijj > 1. Therefore 

V 2 ^ = —KijK* 

is a supersolution, i.e., it satisfies V 2 ^ — ip) < and 
%jj = 1 is a subsolution, i.e., V 2 (-0 — -0) > 0> aiK h of course 
ip > 0. Finally, the solution with linear momentum P = 
lies between. Therefore, as we change P the solution is 
trapped between the sub- and supersolution. The su- 
persolution diverges as r~ 4 as r — ► 0, proving that the 
true solution cannot have any divergence stronger than 
r -4 , and in particular that there are no logarithmic di- 
vergences. This allows us to posit an ansatz for ip consis- 
tent with the allowed blow-up powers, and then check by 
consistency with the full Hamiltonian constraint which 
of those survive to the full solution; and this leads to the 
r~ x / 2 behavior determined in Section HTTl This completes 
our outline of an existence proof, which holds for single 
and multiple-black-hole solutions. The only complication 
arises when the trumpet has angular momentum, but we 
will deal with this case in Section [Vl A more rigorous 
proof along the lines of that for the wormhole-puncture 
case (35l . l37| remains to be constructed, and would be an 
interesting topic for future work. 

Having proved that solutions to this system exist and 
are unique, we now must find them numerically. One po- 
tential problem that is apparent from Eqn. (fT9|) is that 
the Di source terms involve half-integer powers of r near 
the puncture, which affects the accuracy of the elliptic 
solver. Concretely, the Dq term will lead to a r 3 / 2 contri- 
bution to the solution, which we expect to limit the solver 
to 1.5-order accuracy near the puncture, and the D\ term 
will lead to a r 1 / 2 contribution, which we expect will limit 
the solver the 0.5-order accuracy near the puncture [4(| . 
and appears at a lower order than the 7-V2-1/2 ^ r ^ erm 
that we have already accounted for in the ip s i ansatz. 

These expectations are borne out in our results. Fig. [3] 
shows the convergence behavior of the Li norm for the 
entire solution. We find that the convergence is at less 
than first-order, consistent with the half-order conver- 
gence predicted above. (Since we no longer have an an- 
alytic solution to compare with, we evaluate the conver- 
gence by comparisons between solutions with successive 
numbers of collocation points. We chose to sample N in 
multiples of four, and therefore display the Li norm of 
(wjv+4 — ujy) in the figure.) However, if we include in 
the Li norm only that part of the computational domain 
that is outside the apparent horizon of the black hole (lo- 
cated approximately at r = 0.77m), then the errors show 
exponential convergence up to about N = 32. For higher 
numbers of collocation points the convergence rate de- 
teriorates, and for the larger values of N shown in the 
figure the results arc consistent with fourth-order con- 
vergence. This demonstrates that the behavior near the 
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FIG. 3: The error behavior of the Hamiltonian-constraint 
solution for a single boosted trumpet. The dashed line shows 
the convergence of the L2 norm over the entire domain, while 
the solid line shows the L2 norm for the region of the domain 
outside the black-hole horizon. See text for more details. 



puncture limits the accuracy of the solution, but that this 
limitation is essentially localized within the black hole. 

The D\ term is due to the A$/r 5 term in A^A^ 
(which is in turn due to the cross-term between the 
Schwarzschild and Bowen-York extrinsic curvatures). If 
we remove these cross terms from the source function, 
we obtain the convergence behavior shown in Fig. [H 
we now see, as expected, that for N > 32 the conver- 
gence approaches 1.5-order over the entire domain, con- 
sistent with the earlier discussion. Unfortunately, this 
solution does not represent the correct conformal factor 
for a boosted Bowen-York trumpet puncture! 

Although the inclusion of the Bowen-York extrinsic 
curvature limits the accuracy of our solver near the punc- 
ture, the solution is still very accurate over most of the 
computational domain, and is anyway accurate enough 
for most practical purposes everywhere. If one wished to 
produce yet more accurate solutions, one option would 
be to use a coordinate transformation from r to X that 
lead to the solution near the puncture being expanded in 
powers of r 1 / 2 . However, for the purposes of this paper, 
such accuracy is not required, and we simply make this 
observation for future use. 



V. SINGLE SPINNING BOWEN-YORK 
TRUMPETS 

The construction of a solution for a single spinning 
Bowen-York puncture trumpet is complicated by the fact 
that the Bowen-York extrinsic curvature for a spinning 
black hole diverges as 1/r 3 near the puncture. In this 
case, the behavior of the conformal factor near the punc- 
ture will be modified by the presence of the Bowen-York 
term. However, we will show that it is possible to de- 
termine the angular dependence of the divergent term in 
the conformal factor from a local one-dimensional ODE, 
which can be easily solved to construct the appropriate 
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FIG. 4: The error behavior of the Hamiltonian-constraint 
solution for a single boosted trumpet, with the r term re- 
moved from the source term. The solution displays clean ex- 
ponential convergence up to about TV = 32, and then the 
convergence deteriorates to 1.5-order (see text). 



ansatz for a full numerical solution. 



A. Angular dependence at the puncture 

For convenience in what follows, we will express the 
problem in spherical coordinates. We will assume that 
the conformal factor now behaves as ifj ~ D(B)/y/r as 
r — > 0. The square of the conformal extrinsic curvature 
that appears in the Hamiltonian constraint is now 



A 2 = A ij Ai 



6C7 2 

r 6 



18S 2 (l-cos 2 



(20) 



where S is the angular momentum of the black hole. Note 
that in this case there are no cross terms. 

To extend our earlier existence proof to the spinning 
case, we need to take into account the change in the di- 
vergent term in the conformal factor. We start by prov- 
ing a monotonicity condition for the Hamiltonian con- 
straint with these boundary conditions. More precisely, 
Let A\ = A 2 {C, Si) and A\ = A 2 {C, S 2 ), where S 2 > Si. 
Since A\ > A 2 , then tp2 This means that, if we fix 

C and pump up S, the conformal factor monotonically 
increases. 

The proof is as follows. Subtract the two solutions to 
get 

V 2 (0 2 - Vi) + \[A 2 ^ 2 7 - A 2 ^ 7 ] = 0. 

Now multiply across by r m where to lies between 1/2 and 
1, and find an equation for £ = r m {ip2 ~ ^l)- We find 
that 



~ i} to „ m — m 

V 2 £ cU — 

r r 



£+?—{A 2 i> 2 7 -A 2 ^ 7 }=0. (21) 



We can see that £ vanishes both at r = and at infin- 
ity The quantity £ can never be negative because, if it 
were, it would have a negative minimum, and this cannot 



happen. Let us assume that it docs have such a nega- 
tive minimum. Let us see what happens to Eqn. (|2ip 
at that point. We have V 2 £ > 0, — m/r<9 r £ = 0, — (to — 
to 2 )/V 2 £ > 0, and +r m /%[A 2 iq 7 - Afipf 7 ] > 0. The 
last term is the only slightly tricky term. If £ < 0, then 
tp2 < "01 an d 4>2 7 > iPi 7 ■ Since we assume A\ > A\, this 
term is also non-negative and the sum cannot add up to 
zero. 

Now we want to consider how D{9) behaves, where we 
assume ijj = D(0) / \[r + 0{*Jr) near the origin. When we 
substitute into the Hamiltonian constraint, we get the 
following equation for D{&): 



D' 



1 



1 



tan(0) 



8D 7 



D" + — --D + ^ [6C* 2 + 18,S 2 (1 - cos 2 6)] = 0, 

(22) 

where D" is second derivative with respect to 9. This is 
defined on the interval < 9 < tt, but will be symmetric 
around ir/2. At a maximum we have 

D 8 < 1/2[6C7 2 + 185 2 (1 - cos 2 6)}, 

while at a minimum we have 

D 8 > 1/2[6C 2 + 185 2 (1 - cos 2 6)]. 

Therefore the maximum should occur at 7r/2 and the 
minimum at 9 = and D satisfies 

3C 2 < D 8 < [3C 2 + 95 2 ]. 

These upper and lower bounds allow our earlier existence 
proof to go through unchanged. 

Eqn. ([22|) should be read as a one-dimensional second- 
order equation for D{6) on the interval < 9 < ir/2, 
with Neumann boundary conditions, i.e., D' = at both 
ends. 

A solution of Eqn. ([22]) provides the necessary informa- 
tion to construct a single spinning Bowcn-York puncture 
trumpet. The most important feature of Eqn. (|22[) is that 
it is local: we need only solve a simple one-dimensional 
ODE in order to calculate the requisite boundary infor- 
mation — regardless of the linear momentum of the black 
hole, and regardless of the presence or otherwise of other 
black holes in the data. 



Solution of the nonlinear angular-dependence 
ODE 



We solve the nonlinear ODE Eqn. (|22|) by linearizing 
and solving iteratively. The average of the upper and 
lower bounds is used as an initial guess. A simple appli- 
cation of the NDSolve function in Mathematica suffices to 
produce an accurate solution. The solution for S/M 2 = 1 
is shown in Fig. [51 the function D(9) is seen to lie well 
within the upper and lower bounds derived in the previ- 
ous section. 

Figure [5] shows the maximum value of D(9), which oc- 
curs at 9 = 7r/2, as a function of the angular momentum 
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FIG. 5: The solution D(8) for C 2 = 27/16 and S = 1.0. 
The upper and lower bounds, Df omer = 3C 2 and D^ pper = 
3C 2 + 8S 2 are shown with dashed lines. 
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FIG. 6: Values of the maximum value of D(9) (at 9 — tv/2) 
as a function of the angular momentum S, shown with a solid 
line. Also shown as a dashed line is the upper bound. The 
maximum behaves as S 1 ^ 4 for large S. 



S. The figure shows the upper bound on the solution, 
(3C 2 + 9S* 2 ) 1 / 8 , for comparison. The maximum behaves 
as expected, i.e., grows as S 1 ^ 4 for large S. When S is 
small, the 3C 2 term dominates, and the value approaches 
the Schwarzschild value of a/3/2. 

Now that we have calculated D(9), we are able to solve 
the Hamiltonian constraint for both boosted and spin- 
ning Bowcn-York trumpets. The conformal- factor ansatz 



is now provided by replacing the y // 3M/2r term in (|13|) 
with D{6)/tJt. In the numerical procedure to solve the 
Hamiltonian constraint, the derivatives of D(8) required 
in the construction of V 2 ?/^ arc trivial to calculate in 
our Mathematica-based solver, because D(9) is available 
from the solution to (|22p as an InterpolatingFunction 
to whatever precision is required. 

Note, however, that for the spinning case we do not 
know the next-to-leading order behavior of the solution to 
the Hamiltonian constraint (the coefficient of the j-v2-i/ 2 
term) as we did in the boosted case, and this will restrict 
the accuracy of our solver to that given by the ipsi ansatz 
in Section IIVI and of course the magnitude of this term 
will grow with the value of the angular momentum. For 
this reason, high accuracy is difficult to achieve for ex- 
tremely high values of the spin. For the data sets stud- 
ied in this paper we consider angular momenta no higher 
than 5" = 10M 2 , which corresponds to S/M 2 w 0.924. 
We will now discuss the junk radiation content of our 
data sets in more detail. 



VI. RADIATION CONTENT OF 
TRUMPET-PUNCTURE DATA 

Bowen-York black holes can be considered as Kerr or 
boosted Schwarzschild black holes, plus some unphysical 
radiation content, which either falls into the black hole 
or radiates away as junk radiation. We can estimate the 
radiation content of the data as [13, HH 

Erad = \]e\ dm -P 2 - M. (23) 

To evaluate this quantity we first need an estimate of 
the black hole's mass M. The standard way to calculate 
this is via the area of the apparent horizon of the black 
hole. We calculate the irreducible mass, Mi rr = y/A/16n 
and then use the Christodoulou formula [4!| to estimate 
the total mass of a black hole with angular momentum 
S, 

M 2 = Ml r + -^. (24) 

irr 

For boosted wormhole data, the black-hole mass can 
also be estimated by calculating the ADM mass at the 
extra asymptotically flat end; one can see by performing 
an inversion transformation on the Bowen-York extrinsic 
curvature that its contribution at the extra end falls off 
as ?'~ 4 , and therefore we expect that it contributes very 
little junk radiation in the second copy of the exterior 
space. This suggests that the ADM mass evaluated at 
the second asymptotically flat end (i.e., at the puncture) 
will provide a good measure of the mass, and this has 
been confirmed by numerical observations [HI, and the 
"ADM puncture mass" has become a standard tool in 
wormhole puncture data [3, l3~i liol. Iill| . 

There are two drawbacks of the ADM puncture mass. 
One is that it does not provide a good estimate of the 
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mass for spinning black holes, since in that case the 
Bowen-York extrinsic curvature has the same fall-off be- 
havior at both asymptotically flat ends, 0(r~ 3 ), and con- 
tributes roughly the same junk radiation into both ex- 
terior regions. We have verified this in numerical tests, 
where we find that the ADM puncture mass for spinning- 
Bowen- York-puncture data sets equals the ADM mass 
calculated at spatial infinity to within the numerical ac- 
curacy of the solver (w 10 -8 ). 

The other disadvantage of the ADM puncture mass, 
which applies in general to wormhole puncture data, is 
that the mass cannot be prescribed a priori, because 
the relationship between the mass parameter m and the 
black-hole mass M is nonlinear. In order to construct 
Bowen-York wormhole punctures with specific masses, 
an iteration procedure must be used. 

The situation appears to be quite different in the trum- 
pet case. Here the mass parameter m does seem to pre- 
scribe the mass of the black hole, at least for boosted 
black holes. This is presumably related to the fact that 
the Bowen-York extrinsic curvature does not affect the 
geometry of the trumpet, irrespective of the value of the 
linear momentum. This interesting (and useful) prop- 
erty of the boosted Bowen-York trumpet deserves further 
study. 

The same cannot be the case for spinning black holes, 
however, where the coefficient of the singular term in the 
conformal factor is an angular function of the spin. We 
could propose a mass based on the area of the trumpet, 
but this is not necessarily useful, because we do not know 
the relationship between the trumpet area and the black- 
hole mass for spinning black holes. For spinning black 
holes we must make use of the mass calculated from the 
area of the apparent horizon, Eqn. (|24D . 

We are now in a position to estimate the junk radiation 
content of our boosted and spinning trumpet data sets. 

Fig. [7] shows the estimate of the radiation content for 
boosted wormhole and trumpet initial-data sets. We see 
that the results are almost identical for both classes of 
initial data. This also provides further evidence of the 
equivalence of the mass estimates that were used for each 
class of data. These results can further be compared with 
those for other families of boosted Bowen-York data 
53,51], f° r wmcrl the values of the junk radiation content 
appear to be very similar. 

Fig. [5] shows the same quantity estimated for spin ning 
trumpet data sets. If we compare with the results in [43, 
53 we see that the use of the trumpet topology does not 
noticeably change the junk radiation content. 

It was pointed out in [26l . l50l ] that taking the limit as 
in — > while keeping S fixed is equivalent to keeping m 
fixed and taking the limit S — ► oo. In other words, by 
simply removing the Schwarzschild trumpet term from 
the extrinsic curvature, we can construct data equivalent 
to the S — > oo limit. Furthermore, since we know that 
the horizon is located at the puncture for these data, we 
can directly calculate the apparent-horizon area to high 
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FIG. 7: Estimate of the radiation energy content of boosted 
black-hole initial-data sets. The grey squares indicate worm- 
hole data, and the black circles indicate trumpet data. The 
results for both Bowen-York trumpets and wormholes are 
shown. The results are identical at the level of accuracy of 
the data: as one might expect, the use of a trumpet versus 
a wormhole topology does not affect the radiation content of 
the data. 



accuracy from our angular function D{6): 

A = lim / V>V sin(0)d0# (25) 

r->0j 

= 2ir [ D i {e)sin(9)d9. (26) 



We do this and find that SIM 2 — 0.9837, in precise agree- 
ment with the results in [27| . although we note that via 
Eqn. (f2"2")l one can calculate this value to arbitrary ac- 
curacy. We also find that S/M 2 A ri M = 0.928, again in 
agreement with the results in [27] . These numbers pro- 
vide upper and lower bounds on the spin of the final Kerr 
black hole, after the junk radiation has left the space- 
time. We evolved these data, and found that less than 
0.05 % of the energy in the initial slice was radiated away, 
and therefore the rest of the junk radiation falls into the 
black hole (in agreement with the observations in [2(|), 
and the final Kerr black hole has a spin parameter of 
0.928 < a/m < 0.929. Note also that it follows from the 



results in 2(| |50| that the high-angular-momentum lim- 
its of the wormhole and trumpet Bowen-York data are 
equivalent. 



VII. BINARY TRUMPETS 

We now wish to construct data for two Bowen-York 
trumpets. The linearity of the momentum constraint 
with K = allows us to superimpose any number of 
solutions: for each black hole we simply include both 
the Schwarzschild trumpet extrinsic curvature and the 
Bowen-York extrinsic curvature to obtain a valid solution 
of the momentum constraint. For black holes located at 
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FIG. 8: Estimate of the radiation energy content of spin- 
ning black-hole-trumpet initial-data sets, including the ex- 
treme limit, at which S/M 2 = 0.9837 and E rad = 0.0296, i.e., 
the junk radiation never consists of more than 2.96% of the 
energy of the spacetime. 



Y\ and r 2 , the extrinsic curvature is therefore 



A, 



i4£(r-n)+4g*(r-r i; P 
+M(T-r 2 ) 



iff (r-r 2 ;P 2 ) 



(27) 



We once again need a suitable ansatz for the conformal 
factor. The first obvious choice is to generalize the ansatz 
used for a single black hole and try 



ip 



guess 



ioi(n) 



Roi 
n 



R 



3/2 
01 



wi{r 2 ) 



' Ro2 

T 2 

-w 2 (ri)w 2 (r2), 



3/2 



mi 



m 2 



Roi 

7f 

R02 

7f 



(28) 



where i?oi = 3mj/2, p = — 1 — \/2 and q = \/~2 — 1/2, and 
where the w 2 weightings are multiplied so that the result- 
ing function is zero at each puncture, and asymptotes to 
unity far from the source. 

We saw in Sec. IIVI that the requirement that u = 
at the puncture removed any logarithmically divergent 
terms from the solution, but this was possible only be- 
cause the problematic part of the source term was linear 
in u: setting u = removed that term. In the binary 
case, with the ansatz we have chosen, this is not neces- 
sarily so simple. Near one puncture (let us choose r\) the 
conformal factor behaves as ip — ^3m/2r + F + u, where 
F is the contribution from the second term in Eqn. (|28[) . 
This could also generate a logarithmic term. One solu- 
tion would be to determine the appropriate value of u at 
the puncture such that this term no longer contributes 
(i.e., u(r\ = 0) = —F), and enforce this in the solver, 
or hope that the solver finds that value 1 . An alterna- 



1 Research performed concurrently with that in this paper found 
that indeed the solver does appear to locate this value l5ll 




FIG. 9: The function u for the binary configuration described 
in the text, represented in the bi-polar coordinates used in the 
pseudospectral solver. Note that X — 1 corresponds to spatial 
infinity, while X = — 1, Y = ±1 are the puncture locations. 



tive solution is to choose an additional weighting factor 
so that in fact A = 0, and to again impose our standard 
u = or u' = boundary condition at the puncture; this 
is the approach that we will follow. 

To solve the Hamiltonian constraint numerically for 
binary trumpets, we adopt similar coordinates in our 
pseudospectral solver as developed in [23| for use with 
wormhole puncture data. For an equal-mass binary with 
punctures located on the x-axis at x = ±6, we make the 
coordinate transformation 



2b (5 + X(2 + X))Y 
(1 + Y 2 ){3 + X){X - 1) ' 
46(1 + X)(Y 2 - 1) cos <p 
(l+Y 2 ){X 2 + 2X -3) ' 
Ab{l + X){Y 2 - l)sin0 
(l+y 2 )(X 2 + 2X-3) ' 



(29) 
(30) 
(31) 



In these coordinates X = 1 corresponds to spatial in- 
finity. The points (X, Y) = (— 1,±1) correspond to the 
puncture locations at x = ±6. The line along the a;- 
axis between the two punctures is mapped to the plane 
X = — 1. For a full description of this coordinate system 
and its properties, the reader is referred to [23j |. 

These coordinates make it particularly simple to ap- 
ply additional weighting factors that remove at each 
puncture the contribution to the conformal factor ansatz 
from the other puncture. The weights we choose are 
u'i wi cos[(7r/4)(l±F)] 4 . 

As an example, we construct data for the same config- 
uration as in the "D10" case studied in [52} : the punc- 
tures are located at x = ±5M, and the momenta are 
P = (T9-80376 x 10" 4 , ±0.0961073, 0). The specific mo- 
menta are not important for this test; we simply choose 
the same numbers to allow a direct comparison of the 
initial-data sets. 

The solution u for this system is shown in Fig. [SI rep- 
resented in the coordinates ([29)1 - ([3"Tj) . along the plane 
z = (0 = 0). 

The convergence of the solver for these data is shown 
in Fig. 1101 The results indicate surprisingly good conver- 
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FIG. 10: Error behavior for a binary configuration. The 
plot shows the L2 norm of the difference between solutions 
produced with N and N + 4 collocation points. Only values 
along the 2 = plane are included in the calculation, but since 
the punctures lie in this plane, this plot shows the dominant 
error behavior for the solution. 

gence in comparison to the singlc-black-hole cases. This 
may be due to a cancellation in some other problematic 
terms in the binary case. For example, far from the bi- 
nary the Hamiltonian constraint source term will closely 
resemble that of a single spinning black hole; similar can- 
cellation affects may play a role throughout the compu- 
tational domain. 

We can calculate the ADM mass of the system by 
noting that near spatial infinity, ip ~ 1 + Madm/^t), 
and obtain the ADM mass from the radial derivative 
of u as r — > 00. As an indication of the accuracy of 
our solver, and of the level of difference between worm- 
hole and puncture binary data, the ADM masses for the 
wormhole and puncture data with the same choice of 
black-hole mass, separation and linear momenta, were 
0.9897136 and 0.989706, respectively. If we calculate the 
binding energies (E^ = Madm — Mi — M2) for these 
two data sets, they are i?&,wormhoie = —0.0102864 and 
Eb, trumpet = —0.0102939. This demonstrates that these 
data sets are physically extremely close — with the added 
advantage in the trumpet case that the black-hole masses 
could be specified directly through the mass parameter, 
while in the wormhole case they had to be calculated by 
a nonlinear iteration procedure [53 . l52| . 

VIII. NUMERICAL EVOLUTION OF THE 
DATA 

Having proposed and produced a new class of black- 
hole initial data, and claimed certain gauge and physi- 
cal properties for them, we now need to evolve a set of 
trumpet binary data and put our claims to the test. In 
particular, there are two questions we wish to answer: 

1. We expect that the trumpet data are in coordinates 
closer to those preferred by the moving-puncture 
method than wormhole data; is this true? 



2. Do the wormhole and trumpet data describe the 
same physical situation, or, in practical terms, do 
they produce the same gravitational-wave signal? 

We evolve the data using the same version of the BAM 
code [13, HH used to produce the results in [52|, with 
which we compare the gravitational waveform. In the no- 
tation of those works, we use the same TV = 64 grid lay- 
out as used for the "D10" simulation; see Table 1 in [52l |. 

A. Gauge changes 

The first question that we have posed above is difficult 
to answer. The data that we have produced are maxi- 
mally sliced, while in the moving-puncture method one 
usually deals with 1+log slicing, and the data will quickly 
cease to be maximally sliced and will asymptote to their 
appropriate 1+log form. In addition, the punctures are 
initially stationary, but will pick up speed once the evolu- 
tion begins; this constitutes yet another change of gauge. 
These gauge changes may be "larger" than those induced 
by the transition of wormhole data to trumpet form — 
whatever "larger" means in the context of gauge changes. 

However, we can perform one simple test to quantify 
the change in gauge between the two sets of data. In 
wormhole data, the apparent horizons of the two black 
holes are located on surfaces with coordinate radii close 
to r ~ to/2, where to is the mass parameter in the worm- 
hole puncture conformal ansatz (JSJ). For trumpet data, 
on the other hand, the horizon is at about r w 0.78A/. 
If we evolve both wormhole and trumpet data with a 
variant of 1+log slicing that will asymptote to maximal 
slicing for a stationary spacetime, then we expect that 
the horizon radius will stay roughly fixed in the trumpet 
case, while in the wormhole case it will increase quickly 
to a value close to r w 0.78M. (The rapid expansion of 
the horizon early in simulations is standard in moving- 
puncture simulations; sec, for example, @,H|3|.) 

The slicing condition that approaches maximal slicing 
for a stationary solution is 

d t a = -2aK, (32) 

i.e., the standard 1+log slicing used in moving-puncture 
simulations, but without the shift term on the left-hand 
side. With this gauge condition the data will deviate 
from maximal slicing at early times, but will again be 
approximately maximally sliced after about t = 10 M 
of evolution [l{|. In addition we set r\ = in the 
T-driver shift condition, to minimize additional gauge- 
related growth in the horizon [ijj [34[. The results are 
shown in Fig. 1 1 H and are as expected: in the wormhole 
case the horizon radius grows to about 0.75M within 
10M of evolution, while in the trumpet case the horizon 
radius remains close to that value at all times. The addi- 
tional oscillations may be due to other gauge effects, but 
are of much smaller magnitude than the main effect we 
have just described. 
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FIG. 11: Coordinate radius of the apparent horizon as a func- 
tion of time, for one of the black holes in a binary evolution. 
The data are initially maximally sliced. At early times the 
slicing will deviate from K = 0, but in a stationary situa- 
tion would return to maximal slicing within t fa 10M. As 
expected, the apparent-horizon radius shows much less devi- 
ation for trumpet data (solid line) than for wormhole data 
(dashed line). 



B. Junk radiation 
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FIG. 12: Junk radiation pulse from binary wormhole and 
trumpet data, with initial separation of D — 10M, and 
radiation-extraction radius of R^x = 90M. The trumpet- 
data results are shown with the thick dashed line, and the 
wormhole-data results are shown with the grey continuous 
line. As suggested by Fig. [7] the junk-radiation content of 
both data sets is of comparable magnitude. 



property that we claim for a new initial-data set actually 
holds. 



We construct two sets of binary initial data (wormhole 
and trumpet) for a binary with initial coordinate sepa- 
ration of D = 10M. We adjust the initial momenta such 
that both sets of data exhibit quasi-circular inspiral. (We 
find that different values of the initial momenta are re- 
quired for each class of data; the reasons for these small 
differences are at least partially due to the coordinate 
change made manifest by the different apparent horizon 
sizes mentioned previously, and deserve further investi- 
gation in future work.) We then evolve using standard 
moving-puncture gauge choices, i.e., the full 1+log slic- 
ing condition, (dt — f3 1 di)a = —2aK, and with r)/M = 2 
in the T-driver condition. We now wish to evaluate the 
differences in the gravitational-wave signal between sim- 
ulations using each data set. 

The first point of comparison is the burst of junk ra- 
diation at the beginning of the simulation. Based on the 
results in Section HVl we would expect that the junk radi- 
ation is the same in wormhole and trumpet data. Fig. 1121 
shows the pulse of junk radiation in the spin-weight -2, 
(£ = 2, m = 2) mode of rv& 4 , as calculated R ex = 90M 
from the source. (Full details of the wave-extraction pro- 
cedure used in the code are given in [34).) Although the 
junk pulses from the two data sets are not identical, they 
are very similar; it is certainly not possible to definitively 
claim that one type of data contains less junk radiation 
than the other. 

We emphasize that this result is not merely a demon- 
stration of a result that we know to be true. The estimate 
of the radiation content of the initial-data sets, based on 
the initial data alone, is no more than that: an estimate. 
It is only by evolving the data in a full general-rclativistic 
simulation that we can be certain that this (or any other) 



C. Inspiral-merger-ringdown signal 

We now consider the full inspiral-merger-ringdown GW 
signal generated by the inspiral and coalescence of the 
two black holes. In this simulation the binary completes 
about five orbits before merger. 

We focus of the dominant (£ = 2, m = 2) spin- 
weighted spherical harmonic mode of r\E , 4, as extracted 
at R ex = 90M from the source. Fig. [T3] shows sepa- 
rately the inspiral and merger-ringdown portions of the 
real part of H&4 22- (The plot begins after the junk ra- 
diation has passed through the R ex = 90M radiation 
extraction sphere.) The time has been shifted so that 
the maximum amplitude occurs at t = 0. 

The figure includes both the wormhole- and trumpet- 
data results. The results are indistinguishable, except 
for a very small amount of de-phasing early in the sig- 
nal, due to the slightly different effective choice of initial 
parameters. 

Fig. [Til shows the amplitude of r^f ^22 for the wormhole 
and trumpet data. In this case the lines can be distin- 
guished due to the slightly different eccentricities present 
in the two data sets. Once again it is clear, however, that 
the two waveforms agree extremely well; they certainly 
agree well within the error levels discussed in the recent 
Samurai project [H, which demonstrated that waveforms 
that agree to this level are well within the accuracy re- 
quirements for detection and parameter estimation with 
first- and second-generation ground-based GW detectors. 
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FIG. 13: Comparison of the inspiral and merger waveforms. 
The trumpet-data results are shown with the thick dashed 
line, and the wormhole-data results are shown with the grey 
continuous line. A time and phase shift have been applied so 
that the amplitude maxima occur at the same time, at which 
time the waveforms are in phase 



IX. DISCUSSION 

In this work we have extended the puncture method to 
produce trumpet data for boosted, spinning and binary 
black holes based on the Bowen-York extrinsic curvature. 
In the boosted case the generalization is straightforward, 
and in the spinning case a simple one-dimensional nonlin- 
ear ordinary differential equation must be solved to deter- 
mine the angular dependence of the asymptotic trumpet 
geometry. 

We have discovered one surprising advantage of trum- 
pet data over their wormhole counterpart, which is that 
the mass of a boosted Bowen-York trumpet can be 
prescribed analytically by the mass parameter in the 
conformal-factor ansatz used to solve the Hamiltonian 
constraint. This is a great computational advantage over 
the wormhole case, where the mass parameter must be 
iterated to produce data that contain black holes with 
specific desired masses. This relationship could not how- 
ever be extended to spinning black holes. 

The motivation to produce black-hole initial data in 
trumpet form is that this is the topology that is pre- 
ferred by the gauge conditions that are used in the 
moving-puncture method, which is itself the most popu- 
lar method for simulating black- hole binaries. Although 
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FIG. 14: The amplitude of ^4,22, as calculated from moving- 
puncture simulations of wormhole and trumpet puncture ini- 
tial data, with initial separation of D — 10M. The waveforms 
are shifted such that the maxima in the amplitude occur at 
the same time. The thick dashed line shows the trumpet-data 
results, while the continuous grey line shows the wormhole- 
data results. 



we do not expect (and did not find) any dramatic dif- 
ferences in the properties of black-hole simulations be- 
tween wormhole and trumpet data, the construction of 
these data are an important first step towards ideal ini- 
tial data for puncture simulations. Such ideal data will 
be in the 1+log gauge (or whatever slicing condition is 
ultimately used to evolve the data, one natural alterna- 
tive being hyperboloidal slicing conditions [3, HH ) , will 
represent true boosted Schwarzschild or Kerr black holes, 
and will be in trumpet form. In this work we have made 
the simplest step in this direction, i.e., we have produced 
trumpet data, but they are maximally sliced and repre- 
sent only approximations to boosted Schwarzschild and 
Kerr black holes. 

Efforts in these other directions have already been 
made. Data for superposed Kerr punctures have for ex- 
ample been presented in [31(, and superposed boosted 

55j ; non- 



Schwarzschild punctures have been used in 
conformally-flat data that attempt to include the GW 
signal from the earlier inspiral of the binary have been 
proposed in [56|. Work has also been done in produ cing 
non-conformally-flat data with excision techniques [27I . 
\El\ . It is also now known how to produce 1+log trumpet 
puncture data for a single Schwarzschild black hole [F3 |. 
It is likely that a combination of all of these approaches 
will be necessary to produce the optimal data for punc- 
ture simulations. 
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